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Used in the fixed-target mode, the multi-TeV LHC proton and lead beams allow for studies of heavy-flavour 
hadroproduction with unprecedented precision at backward rapidities -far negative Feyman-x- using conven¬ 
tional detection techniques. At the nominal LHC energies, quarkonia can be studied in detail in p ■¥ p, p -v d 
and p -v A collisions at sfstm — 115 GeV as well as in Pb -(- p and Pb -l- A collisions at sfsiff ^ 72 GeV with 
luminosities roughly equivalent to that of the collider mode, i.e. up to 20 fb^'yr^' in p-v p and p-vd collisions, 
up to 0.6 fb^'yr^' in p -l- A collisions and up to 10 nb^'yr^' in Pb + A collisions. In this paper, we assess the 
feasibility of such studies by performing fast simulations using the performance of a LHCb-like detector. 


I. INTRODUCTION 

Since its start-up, the large hadron collider (LHC) -the most 
energetic hadron collider ever built so far- has already made 
the demonstration of its outstanding capabilities. These can 
greatly be complemented by the addition of a fixed-target 
physics program. Its multi-TeV beams indeed allow one to 
study p -\- p, p -\- d and p -\- A collisions at a center-of-mass 
energy yfsNN - 115 GeV as well as Pb -i- p and Pb -i- A col¬ 
lisions at - 72 GeV, with the high precision typical of 
the fixed-target mode. In this context the proposal of a fixed 
target experiment at the LHC [1], referred to AFTER®LHC - 
A Fixed Target Experiment -, has been promoted [ 1 ] in order 
to complement the existing collider experiments such as the 
Relativistic Heavy Ion collider (RHIC) or the future Electron- 
Ion Collider (EIC) project in a similar energy range. The idea 
underlying the AETER@LHC proposal is a multi-purpose de¬ 
tector allowing for the study of a multitude of probes. 

Various technological ways to perform fixed-target exper¬ 
iment at the LHC exist. On the one hand, the beam can be 
extracted by means of a bent crystal. This technology [2, 3] 
is currently developed as a smart beam-collimation solution 
and is studied by the UA9/LUA9 collaboration respectively at 
SPS and LHC. A bent crystal installed in the halo of the LHC 
beam would deflect the particles of the halo onto a target, with 
a flux of 5 xlO*^ proton/s without any impact on the LHC per¬ 
formances [3-5]. 

On the other hand, the LHC beam can go through an 
internal-gas-target system in an existing (or new) LHC exper¬ 
iment. Such a system is already tested at low gas pressure by 
the LHCb collaboration to monitor the luminosity of the beam 
[6-8]. Data were taken at a center-of-mass energy of s[snn - 
87 (54) GeV with p -H Ne (Pb -H Ne) collisions during pilot runs 
in 2012 and 2013. Although this system, called SMOG, was 
tested during only few hours in a row, no decrease of the LHC 
performances was observed. 

In the bent-crystal case, the luminosity achievable with AL¬ 
TER @LHC would surpass that of RHIC by 3 order of magni¬ 


tudes [1]. We have reported in Table (I) the instantaneous and 
yearly integrated luminosities expected with the proton and Pb 
beams on various target species of various thicknesses, for the 
bent-crystal as well as internal-gas-target options. Integrated 
luminosities as large as 20 fb ' can be delivered during a one- 
year run of p -I- H collisions with a bent crystal. Besides, it 
is worth mentioning that both technologies allow one to po¬ 
larise the target, which is an important requirement to lead an 
extensive spin physics programme [1,9]. 

Overall, thanks to the large luminosity expected, AL¬ 
TER® LHC would become a quarkonium [10], prompt pho¬ 
ton and heavy-flavour observatory [1, 1 1] in p -n p and p H- A 
collisions where, by instrumenting the target-rapidity region, 
gluon and heavy-quark distributions of the proton, the neu¬ 
tron and the nuclei can be accessed at large x and even at 
X larger than unity in the nuclear case [12]. In addition, 
the fixed-target mode allows for single-target-spin-asymmetry 
measurements over the full backward rapidity domain up to 
Xp ^ -1 [13, 14]. Also, the versatility in the target choices of¬ 
fer a unique opportunity to study the nuclear matter versus the 
hot and dense matter formed in heavy-ion collisions which 
can be studied during the one-month lead run. In the latter 
case, modern detection technology (such as high granularity 
calorimeter) should allow for extensive studies of quarkonium 
excited states, from the tf/{2S) to the Xc and Xb resonances 
thanks to the boost of the fixed-target mode [15]. 

In this paper, we report on a feasibility study of quarkonium 
production at a fixed-target experiment using LHC beams. In 
section II, we outline the simulation framework which was 
used. In section III, we describe how a fast simulation of a de¬ 
tector response has been implemented, following a LHCb-like 
detector setup. In section IV, we present the charmonium and 
bottomonium family studies performed with the p H- H simu¬ 
lations at V* = 115 GeV. In section V, we present multiplicity 
studies in p -H A and A H- p collisions as well as the expected 
nuclear modification factors for J/i/r and T in p H- Pb collisions 
at sfstfiv = 115 GeV. Einally in section VI some prospects for 
Pb -H A measurement at yfsNN - 72 GeV are given. Section 
VII gathers our conclusions. 
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TABLE I. Expected luminosities obtained for a 7 (2.76) TeV proton 
(Pb) beam extracted by means of a bent crystal or obtained with an 
internal gas target. 
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II. SIMULATION INPUTS 


In order to get the most realistic minimum bias simula¬ 
tions at AFTER@LHC energy for quarkonium studies in the 
dimuon decay channels, we have simulated the quarkonium 
signal and all the background sources separately to have under 
control the transverse momentum and rapidity input distribu¬ 
tions as well as the normalisation of the different sources. 

The simulation has been performed for p + p collisions at 
V* = 115 GeV. On the one hand, the quarkonium signal and 
the correlated background (Drell-Yan, cc, bb) were simulated 
with HELAC-Onia [16] which produces outputs following the 
format of Les Houches Event Eiles [17]. The outputs were 
then processed with Pythia (Pythia 8.185 [18]) to perform the 
hadronisation, the initial/final-state radiations and the decay 
of the resonances. On the other hand, the uncorrelated back¬ 
ground was obtained from minimum bias p+p collisions gen¬ 
erated with Pythia. 

The relative normalisation of the signal and background 
sources was performed according to the production cross sec¬ 
tion of the process (taking into account initial phase space 
cuts, if any). Values of the cross section and the number of 
simulated events Nsim -not to be confused with the expected 
events for a specific luminosity- are reported in Table (II). The 
cross section values are integrated over rapidity and pj. 


A. Signal and correlated background 


I. Quarkonium signal 


J/i//, i/i(2S), T(15), Y(25) and T(35) were simulated in a 
data-driven way. The amplitude of gg —> (3 H- A (where Q 
is the quarkonium) is expressed in an empirical functional 


form [19]: 


l^^gg^Q+xr - 


Kexp(-Kj^) 


when Pt ^(Pt) 

( 1 ) 

when Pt > (pt) 


where K = A^KsfM^ with s the partonic center-of-mass en¬ 
ergy and Mq the mass of the quarkonium Q taken from the 
PDG table [20]. 
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FIG. 1. Some illustrative comparison between fits and the PHENIX 
data [21] for charmonium production (upper panel) and LHCb 
data [22] for bottomonium production (bottom panel). 

The parameters k. A, n and (px) were determined by fit¬ 
ting the differential cross section cficrldpTdy to the exper¬ 
imental data. The dedicated codes used to perform the fit 
and to generate unweighted events for quarkonium produc¬ 
tion have been implemented in HELAC-Onia [16] and we used 
MSTW2008NLO PDF set [23] provided in LHAPDF5 [24] 

and the factorisation scale Pf - + p^. In order to con¬ 

strain the non-trivial energy dependence of quarkonium pro¬ 
duction, we used the differential measurements of charmo¬ 
nium production performed by the PHENIX collaboration at 
RHIC, in p + p collisions at yfs - 200 GeV [21] to predict 
the corresponding yields at s/s = 115 GeV. Given the lack of 
such measurements for T at RHIC, we performed a combined 
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(T,o, (mb) 

^sim 

JM 

1.30 X 10^3 

1.47 X 10® 

<A(2S) 

1.61 X lO-'* 

1.12 X 10® 

T(1S) 

4.30 X 10^^ 

1.46 X 10® 

T(2S) 

1.22 X 10-'' 

1.49 X 10® 

T(3S) 

5.28 X 10-* 

1.48 X 10® 

Drell-Yan (M > 2.5 GeV/c^) 

2.52 X 10-® 

4.3 X 10® 

Drell-Yan (M > 7 GeV/c^) 

1.49 X 10-' 

2.0 X 10® 

CC 

2.29 X 10-‘ 

81.5 X 10® 

bb 

4.86 X 10-* (gg bb) 

32.3 X 10® (gg bb) 


1.49 X lO-'* (qq bb) 

85.7 X 10® (qq bb) 

minimum bias 

26.68 

11.0 X 10* 


TABLE II. Total cross section for different processes in p + p collisions at = 115 GeV and number of simulated events Aj,, 


fit to CDF [25], ATLAS [26], CMS [27] and LHCb [22, 28] 
data on T production. The values of the htted parameters are 
listed in Table (III). For illustration, the comparison between 
fits and the selected experimental data is shown in Fig. 1 . 
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A 
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<A(25) 
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4 

1.12 

Tfli') 

0.707 

0.0837 
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1883 

T(25') 

0.604 0.0563 

205 

856 

T(35) 

0.591 

0.0411 

197 

886 


TABLE III. Fit parameters obtained after a combined fit of 
cfajdpjdy to the PHENIX data [21] for charmonium production 
and to CDF [25], ATLAS [26], CMS [27] and LHCb [22, 28] 
data for bottomonium production. We have fixed n = 2 and 
(Pt) = 4.5 (13.5) GeV/c for charmonium (bottomonium) production. 
The number of fitted data points is also reported. 

In order to increase the statistics of the simulated data sam¬ 
ple, the decay of the quarkonium in Pythia is forced into 
the dimuon decay channel. The simulated yields are then 
weighted by the cross section for this process multiplied by 
the Branching Ratio (BR). 


2. Open charm 

Open-charm production was simulated with the process 
gg cc in HELAC-Onia. In order to avoid the huge theoretical 
uncertainties in the state-of-the-art perturbative calculations, 
open charm yields at ^fs - 115 GeV are also computed in a 
data-driven way following the method described in the previ¬ 
ous section. Similarly, the matrix element of gg cc is deter¬ 
mined using Eq. (1). The parameters are obtained from a fit to 
the px-differential cc cross section measured by the STAR ex¬ 
periment [29] in pH-p collisions yfs - 200 GeV (see Fig. 3). 
We obtained k = 0.437, A = 3.04 and (px) = 2.86 GeV/c 
when n - 2 hy using CTEQ6L1 [30] and by hxing the c 
quark mass to = 1.5 GeV/c^ and the factorisation scale 


to pf = + p^. The;^^^ of the fit is equal to 4.39 with 10 

experimental data points. The tuned result is shown in Fig. 2. 
The evolution of the cross section with the energy down to 
V* = 115 GeV is then given by HELAC-Onia. 
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FIG. 2. A comparison between fit and the STAR data [29] in p + p 
collisions sA = 200 GeV, for cc production. 


After embedding the Les Houches Event Eile into Pyxhia, 
muons from the underlying Pythia event can be produced on 
top of muons from the initial cc pair. The combination of 
those additional muons with a muon from the initial cc pair is 
not included in our definition of open charm correlated back¬ 
ground. We have however checked that this contribution is 
negligible. In order to increase statistics, the D°, D^, 
and were forced to decay into muons and only those de¬ 
cay muons were considered as correlated background, p^p- 
pairs coming from all possible combinations; D^D^, D^D , 
D^D^, and are considered. The 

simulated events are weighted by the production cross section 
times the pair Branching Ratio times the fraction of c quark 
fragmenting to D°, D^, or This fraction is obtained 
from Pythia and found to be 95%. 
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3. Open beauty 

The theoretical uncertainty on open beauty production is 
relatively smaller than the one on open charm production. 
We therefore calculated open-beauty-production yields with 
a Leading Order (LO) matrix element and which was nor¬ 
malised to the Next-To-Leading-Order (NLO) K factor. The 
NLO cross section with the same setup was calculated by 
MadGraph5_aMC@NLO [31]. We used CTEQ6L1 (CTEQ 6 M) 
for the LO (NLO) calculation. The K factor is found to be 
1.83. The renormalisation and factorisation scale is pr - 

Up - yjml + with the mass of the b quark taken as 

nth - 4.5 GeV/c^. We have adopted a similar definition for the 
open beauty correlated background as the one of open charm 
(see previous section). 


4. Drell-Yan 
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Drell-Yan (DY) correlated background was simulated with 
the process qq y* jZ ^ at LO where qq is a pair of 
same flavour light quarks. The LO calculation was done with 
the CTEQ6L1 pdf set and the renormalisation and factorisa¬ 
tion scale was set to pr - pp - px. In order to have enough 
statistics in the J/^/r and tl/(2S) mass window, a phase space 
cut requesting that the invariant mass of the dimuons (M) is 
greater than 2.5 GeV/c^ was applied. For the simulation of 
the DY background under the T family peaks, a phase space 
cut M > 1 GeV/c^ was applied. The DY cross section ob¬ 
tained with HELAC-Onia at ^fs - 38.8 GeV is compared to 
the existing E 866 data at the same energy [32]. A K factor 1.2 
is needed to match the data and therefore it was also applied 
at ■\fs - 115 GeV. Such a K factor is known to approximately 
account for the higher-order QCD corrections. 

B. Uncorrelated background 

The uncorrelated background was obtained from a min¬ 
imum bias Pythia p + p simulation at yfs = 115 GeV 
using the process So£tQCD:nonDifTractive with the 
MRSTMCal.LHgrid LHAPDF (6.1.4) set [34]. By com¬ 
paring our simulation of open charm with a low statistic 
pure minimum bias Pythia simulation, we have checked 
that the contribution of dimuons originating from a muon 
from charm/beauty and a muon from Tt/K is negligible. The 
dominant source of uncorrelated opposite-sign muon pairs 
is the simultaneous semi-muonic decay of uncorrelated n 
and/or K. In order to avoid possible double counting of signal 
and correlated background processes, the following hard 
processes have been switched off from the minimum bias 
simulations: HardQCDihardccbar, HardQCD: hardbbbar, 
WeakSingleBoson:£fbar2gmZ', Cbarmonium:all and 
Bottomonium:all,. 


* in order to avoid Drell-Yan pair production. 


FIG. 3. Drell-Yan cross section as a function of the rapidity in the 
center-of-mass frame obtained with HELAC-Onia -i- Pythia at yfs = 
38.76 GeV and rescaled by a factor 1.2, together with E866 data 
extracted from [33]. The invariant mass range considered is 7.2 < 
M < 8.7 GeV/cL 


III. FAST SIMULATION OF THE RESPONSE OE A 
LHCb-LIKE DETECTOR 

The HELAC-Onia and Pythia generators provide the 
opposite-sign muon pairs from quarkonia decays, correlated 
and uncorrelated backgrounds sources, as defined in the pre¬ 
vious section. In order to account for the detector resolution 
and the particle identification capabilities of a given detector 
and to investigate the feasibility of the quarkonium studies in 
p + p collisions at Ya ^ 115 GeV, the detector response needs 
to be simulated. For this purpose we have chosen a detector 
setup similar to the LHCb detector [35]. The forward detector 
is very well suited as a fixed-target experiment setup as well, 
with a good tracking and particle identification capabilities. 

According to LHCb analysis cuts, muons in our simulations 
are required to have their transverse momentum satisfying 
Pt > 0.7 GeV/c [36] and their pseudo-rapidity in the labora¬ 
tory frame satisfying 2 < 77 < 5. The 77 cut range corresponds 
to the LHCb detector coverage. Since the momentum resolu¬ 
tion reported by LHCb is 6p/p ~ 0.4 (0.6)% for a momentum 
of 3 (100) GeV/c [37] we consider a momentum resolution of 
dpip - 0.5 %. The single p identification efficiency is taken 
to be ep - 98%, which is an average efficiency obtained by 
LHCb for muons coming from J/i/i decays, for p > 3 GeV/c 
and pt > 0.8 GeV/c [37]. These cuts and the abovementioned 
detector response on the muons are applied to simulate the 
quarkonium states and all the background sources. 

In the case of uncorrelated background, as discussed in sec¬ 
tion II, most of the p originate from or decays. 
If a TT or decays to a 77 before 12 m along the z axis, 
the p is rejected by the tracking system and it is not con¬ 
sidered in the simulation. 12 m corresponds to the distance 
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FIG. 4. Misidentification probability of n (left plot) and K (right plot) as muon candidates as a function of momentum, Pmid(^ —> p) and 
Pmid(K —> yu), respectively. 





FIG. 5. Muon pair, identification efficiency as a function of the pair invariant mass (a), transverse momentum (b) and rapidity (c) 

for uncorrelated muon background. The efficiency takes into account the identification efficiency of the prompt muons, and the n and K 
misidentification probability , > p) and Pmid(K —> /r). 


where calorimeters, followed by muon stations, are placed in 
the LHCb detector setup. If the /t is produced beyond 12 m 
or if a njK is misidentified with yu in the muon stations, a 
jijK misidentification probability is applied. The misiden¬ 
tification probabilities depend on the total particle momen¬ 
tum and were reported by the LHCb collaboration in [38]. 
These probabilities are parametrised with the following func¬ 
tions: Pmid(^ —> p){p) - (0.5 -I- 6.63 exp(-0.13p))% and 
Pmid(K p)(p) - (0.5 -H 8.6 exp(-0.11/7))%, and are shown 
in Fig. 4 (a) and (b), for n and K respectively. Based on the 
single yu identification efficiency e^+/-, the dimuon, p^p^, ef¬ 
ficiency is calculated as a product of the single efficiencies: 
f/i+yi- = X f/y- • For muons coming from or decays, 
misidentification probabilities are used: e^+/- = Putwi^ 
p)ip) or e^+/- = Pmid(K p)(p), respectively for n and K, 
and for prompt muons e^+/- - ep - 0.98. 

The pair efficiency is extracted in each kinematic phase- 
space point and is shown as a function of the dimuon invari¬ 
ant mass, transverse momentum and rapidity in Fig. 5. This 
efficiency is used to correct dimuon spectra obtained with the 
uncorrelated background Pythia simulations. 


IV. QUARKONIUM PRODUCTION STUDIES IN pH 
COLLISIONS AT Vi = 115 GEV 

In this section, we show results on the quarkonium produc¬ 
tion studies in the dimuon decay channels, with the dominant 
background sources. Simulations have been performed for 
a 7 TeV proton beam on a hydrogen target {p + p), which 
gives Vi = 115 GeV. We consider an integrated luminosity 
of 10 fb * which is expected to be obtained after half of a 
LHC year with the crystal mode, as described in section I and 
Table (I). 


A. Background studies 

These simulations allow to quantify the background sources 
in the quarkonium studies in the dimuon decay channel, which 
could potentially make the quarkonium signal extraction more 
difficult or even prevent from obtaining a clear signal. In 
particular, this may be critical for the excited states. We 
present here simulations of invariant mass of opposite-sign 
muon pairs, p^p^, from the quarkonia and from the domi- 
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FIG. 6. Dimuon invariant mass distributions for J/i/r and i^(25), left plot, and Y(nS ), right plot, with different background sources. 



8 10 12 
p!^'* (GeV/c) 


2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 

Y 


FIG. 7. pt (a) and y (b) spectra of J/ij/ signal and different background sources, in the J/tf/ mass range. 


nant background sources, in two mass ranges, see Fig. 6. The 
first range corresponds to the J/ifr and \p{2S) invariant mass 
windows and the second one to the mass range of the T(15), 
T(25) and T(35). The invariant mass distributions are in¬ 
tegrated over the whole transverse momentum and rapidity 
ranges. The plots show the simulated quarkonium signals and 
the background, separately from the different sources, and the 
black solid line is a sum of all contributions. The background 
sources correspond to an integrated luminosity of 10 fb *. 

In the J/i/r and ^{2S) invariant mass window, the dominant 
background source is from uncorrelated pairs, mostly 
from and decays. Contributions from Drell-Yan and 
hh continuum are very small. In the case of T(nS) states the 
Drell-Yan contribution is the dominant one. Under the T{nS) 
peak, a contribution from the cc continuum is negligible, and 
is not considered here. 

The significance {sig = 5/ -\/(S + B), where S is the num¬ 
ber of signal counts and B the number of background counts, 
in the invariant mass range Mq ± 3 <Tq) and the signal to back¬ 
ground ratio (S / B) of each quarkonium state are given in the 


following: 

• = 134.6 10^ cr, S/Bji^ = 4.21 

• si8i/y(2S) = 735.2 cr, SjB^igs) - 0.16 

• ■sl.?T(is) = 140.73 cr, 5 /Bt(is) = 1-75 

• sig’t( 2 s) — 45.29 cr, SIBr{2S) — 0.48 

• — 25.75 cr, S /Bt(3S) — 0.28 

for J/i/r, i/r(25), T(15), T(25) and T(35), respectively. 

Transverse momentum and rapidity distributions for the 
quarkonium signals and for each background source were also 
studied. As an example, pj and y distributions in the }/i// 
mass range, 3.063 < M^+^- < 3.129 GeV/c^ (corresponding 
to Mj/^ + 3 (Tj/^) are shown in Fig. 7. It is visible that the 
distributions for the J/i/z and different backgrounds differ. In 
more backward or forward rapidity regions, the signal to back¬ 
ground ratio increases. This can also be seen in Fig. 8, where 
the dimuon invariant mass distributions in J/tf/ and il/{2S ) mass 
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FIG. 8. Dimuon invariant mass distributions, for three rapidity bins, in J/^/^ and iff(2S) mass window. 2 < < 3, 3 < < 4, 

4 < < 5, shown on panels (a), (b) and (c), respectively. 



FIG. 9. Transverse momentum (a) and rapidity (b) distributions for }/if/, if/(2S), T(15), T(2S) and T(35), from the top to the bottom distribution. 


window are shown in three rapidity ranges. In terms of trans¬ 
verse momentum, one can obtain a very clean signal when 
going to higher pj. Above ~ 4 GeV/c, the uncorrelated back¬ 
ground starts to vanish. Since cc, bb and Drell-Yan simula¬ 
tions are LO simulations, the px spectra of these correlated 
background sources are not shown here. 

B. Quarkonium simulations 

We have also studied the px and rapidity coverage reach of 
the quarkonium signals. The transverse momentum distribu¬ 
tions are shown in Fig. 9 (a), for J/i/r, tf/(2S), T(15), T(2S) 
and T(3S), from the top to the bottom distribution. Similarly, 
Fig. 9 (b) shows the rapidity distribution for each quarkonium 
state. With an integrated luminosity of 10 fb ' the quarko¬ 
nium studies can be carried out in a wide rapidity and px 
range. It should be possible to study T(nS) signals up to 
Px - 10 GeV/c, and J/i/r and i/'(25) could be studied even up 
to Px - 15 GeV/c. All the quarkonium states can be measured 
down to Px = 0 GeV/c. 

This study is limited by the rapidity range of 2 < y < 5, 
in the laboratory frame, due to the pseudorapidity cuts on the 


decay p. The red x-axis on the top of Fig. 9 (b) denotes the 
rapidity in the center-of-mass frame. The rapidity shift for a 
7 TeV proton beam on a hxed target is -4.8, i.e. ycM - ^ ^ 
yiab - 4.8. }/ifr and (1/(28) signals can be studied in the whole 
mentioned rapidity range, while the lowest rapidity reach for 
T(n5) is-2.5-3. 


V. QUARKONIUM MEASUREMENTS IN p+A 
COELISIONS AT = 115 GEV AND Pb+H COLLISIONS 
AT = 72 GEV 

A. Multiplicity in proton-nucleus collisions 

In proton-nucleus collisions, the high track multiplicity 
may induce a high detector occupancy and lead to a reduc¬ 
tion of the detector capabilities. Since LHCb has success¬ 
fully measured the J/ip and T production in pH-Pb collisions 
at = 5 TeV [39, 40], one would expect a good capability 
of such detector under similar particle multiplicity environ¬ 
ment. In the following, the charged particle multiplicity has 
been generated with the EPOS generator [41, 42] in different 
conhgurations; p -i- Pb collisions at yAsw = 5 TeV in collider 
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mode, p + Vh collisions at ^fs^ =115 GeV and Pb + H col¬ 
lisions at = 72 GeV in fixed-target mode. The charged- 
particle multiplicity is dominated by the n multiplicity. By 
comparing these three distributions as a function of the pseu¬ 
dorapidity of the particle in the laboratory frame as shown in 
Fig. 10 , one can conclude that the charged particle multiplic¬ 
ity in a fixed-target mode never exceeds the one obtained in 
p + 'Ph collisions at yisw = 5 TeV in the collider mode in the 
full pseudorapidity range; a detector with the LHCb capabili¬ 
ties will be able to run in such conditions. 



FIG. 10. Averaged number of charged particles m p +A collisions as 
a function of the pseudorapidity in the laboratory frame. 


B. Prospects for the measurements of the nuclear modiflcation 

factors for J/i/^ and T in pH-Pb collisions at =115 GeV 

To illustrate the potential offered by AFTER @LHC in p + 
Pb collisions at yfs^ - 115 GeV, we have evaluated, in this 
section, the impact of the nuclear modification of the gluon 
densities in nucleons within large nucleus -generically re¬ 
ferred to as gluon shadowing- and its uncertainty as encoded 
in the nuclear PDF set EPS09. Eor that, we have used the 
probabilistic Glauber Monte-Carlo framework, JIN [43, 44], 
which allows us to encode different mechanisms for the par- 
tonic production and to interface these production processes 
with different cold nuclear matter effects, such as the afore¬ 
mentioned shadowing, in order to get the production cross 
sections for proton-nucleus and nucleus-nucleus collisions. 
JIN also straightforwardly computes any nuclear modification 
factor for minimum bias collisions or in specific centrality 
classes. In the case of proton-nucleus {p + A) collisions, it is 
the ratio of the yield per inelastic collision in p + A collisions 
to the yield in pp collisions at the same energy multiplied by 
the average number of binary collisions in a typical p + p col¬ 
lision, (Vcoii); 


RpA - 


dN, 


pA 




( 2 ) 


In the presence of a net nuclear effect, RpA is defined such that 
it differs from unity. In the simplest case of minimum bias 
collisions, one should have 


AdCTpp ■ 


(3) 


As in [45], we have used the central curve of EPS09 as well 
as four specific extreme curves (minimal/maximal shadowing, 
minimal/maximal EMC effect), which reproduce the envelope 
of the gluon nPDE uncertainty encoded in EPS09 LO [46]. 

In addition to the modification of the partonic densities, 
quarkonium production in p + A collisions can be affected 
by other effects, for instance by the nuclear absorption which 
depends much on the nature of the object traversing the nu¬ 
clear medium. If the meson is already formed, it may be af¬ 
fected more than a smaller pre-resonant pair. To discuss such 
an effect, it is useful to introduce the concept of the formation 
time, tf, based on the Heisenberg uncertainty principle and the 
time -in the rest frame of the meson- to discriminate between 
two S states, for instance the J/t/r and the il/(2S). In fact, one 
finds [45, 47] that such a time is similar for the charmonium 
and bottomonium states and is on the order of 0.3-0.4 fm. Ob¬ 
viously, this time has to be boosted in the frame where the 
nuclear matter sits. Eor tf smaller than the nucleus radius, the 
quarkonium is formed before escaping it. In the fixed-target 
mode with a proton beam and a nuclear target, the boost fac¬ 
tor is simply jiyiat) - cosh(y;ai,). We therefore obtain tf as in 
Table (IV). 


ycMs 

yiab yiyiab) 

rf'"(y) 

fCMS 

yiab jiyiab) 

-2.5 

2.3 

5 

1.75 fm 

-0.5 

4.3 

37 

13 fm 

-1.5 

3.3 

14 

5 fm 

0. 

4.8 

61 

21 fm 

-1.0 

3.8 

22 

8fm 

0.5 

5.3 

100 

35 fm 


TABLE IV. Boost and formation time in the (target) Pb rest frame of 
a} /if/ and a T as a function of its CMS rapidity at ySiviv = 115 GeV. 

One sees that looking at quarkonium production in p H- Pb 
collisions at different backward rapidities allows one to look at 
quarkonia traversing the nuclear matter at very different stages 
of their evolution. This effect could theoretically be studied 
by giving an ad-hoc rapidity dependence to the effective ab¬ 
sorption cross section, This is left for a future study 

since, here, we wish to consider only the nPDE effects and 
the expected statistics. Other effects to be considered are the 
coherent energy loss [48] (expected to grow in the forward re¬ 
gion) and the rescattering by comovers [49] (expected to grow 
with the multiplicity along the J/i/z direction). 

Since we wish to assess the descriminating power of pos¬ 
sible data to be taken with APTER@LHC, we attribute to the 
EPS09 central values statistical uncertainties which directly 
follow from the differential yields repectively expected in p+p 
and p-i-Pb collisions. Eor that we take an integrated luminosity 
of 10 fb“^ for the p + p runs and 100 pb * for the p -H Pb runs, 
in accordance with the luminosities discussed above (see Ta¬ 
ble (I)). As this stage, we do not consider additional systemati¬ 
cal uncertainties. This simplifying assumption could be lifted 
in a more detail study which would also take into account a 
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FIG. 11. Nuclear modification factor for T as a function of (a) ycMS and (b) pi in p + Pb collisions at = 115 GeV. The uncertainties 
attached to the central points are derived from the statistics to he collected with £.p+p = 10 fb“* and -£p+pb = 100 pb“'. 




FIG. 12. Nuclear modification factor for J/i/r as a function of (a) j'cms and (b) pi in p + Pb collisions at =115 GeV. The uncertainties 
derived from the statistics to he collected with £.p+p = 10 fh“* and -£p+pb = 100 pb“' are smaller than the point size. 



FIG. 13. Nuclear modification factor for J/i^ from h as a function 
of vcMS in p + Pb collisions at = 115 GeV. The uncertainties 
derived from the statistics to be collected with X.p+p = 10 fb“* and 
i^p+Pb = 100 pb“' are smaller than the point size. 


possible detector acceptance (and related efficiencies) as done 
in the previous section. In particular, we do not expect that the 
rapidity region for ycMS >1.5 would be easily accessible. 


In Fig. 11, we show the rapidity dependence of Rp+ph for 
T and its pj dependence near y = 0. The million of T to 
be collected per year allows for the measurement of a Rp+ph 
with a much better precision than the gluon nPDF, nearly up 
to jc —> 1. In addition, one notes that the nuclear modification 
factor is certainly measurable up to px - 10 GeV/c. 

In Fig. 12, we also show the rapidity dependence of Rp+ph 
for i/i/z and its px dependence near y = 0. In both cases, the 
luminosity to be taken in a year at AFTER@LHC yields to 
statistical uncertainties which are largely negligible as com¬ 
pared to the nPDF uncertainties - the statistical uncertainties 
are not even visible on Fig. 12. We except this to hold also for 
the 1^(25) although its yields are down by a factor of 100. 

As aforementioned, the nPDFs do not account for all the ex¬ 
pected nuclear matter effects. However, it is clear that combin¬ 
ing the measurements of T, J/if/ and i/'(25) for -3 < ycMS < 0 
(as a LHCb-like detector would do) will allow one to pin down 
the existence of a possible gluon EMC and antishadowing ef¬ 
fect. We also stress that the complications induced by a ra¬ 
pidity dependence of could be avoided by the parallel 

measurement of Rp+pb for non-prompt J/if/ which can only be 
sensitive to the energy loss since the b quark decay (weakly) 
into the J/i/z, way outside the nucleus. Eig. 13 shows that the 
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trend is similar than for T. Measuring the px dependence of 
Rp+Pb for prompt i/tf/ and T should also avoid the sensitivity 
on formation time effects. 


VI. PROSPECTS OF Pb+A MEASUREMENTS AT 
Vi = 72 GEV 

The charged particle multiplicity has been generated with 
the EPOS generator [41, 42] in different conhgurations; Pb + 
Pb at V^W = 5.5 TeV in collider mode, Pb + Ar, Pb + Xe and 
Pb + Pb at Viw = 72 GeV in fixed-target mode. By com¬ 
paring these three distributions in the pseudorapidity of the 
particle in the laboratory frame as shown in Fig. 14, one can 
conclude that the charged particle multiplicity in a fixed-target 
mode never exceeds the one obtained in Pb -H Pb collisions at 
V^NN = 5.5 TeV obtained in a collider mode in the full pseu¬ 
dorapidity range: a detector with the ALICE MET-tMuon de¬ 
tector [50] capability will be able to run in such conditions. 
Details studies are needed to evaluate up to which multiplicity 
a detector such as LHCb would be able to take good quality 
data. 



FIG. 14. Averaged number of charged particles in A -l-A collisions as 
a function of the pseudorapidity in the laboratory frame. 


VII. CONCLUSION 

In summary, we have shown that in a fixed target mode with 
an integrated luminosity of 10 fb ', using 7 TeV LHC proton 
beam on a hydrogen target, and with a detector setup and per¬ 
formances similar to the LHCb detector, quarkonium studies 
in the dimuon decay channel can be performed over a wide 
transverse momentum range and rapidity in the center of mass 
from-2.8 for Jjij/ and il/{2S), and-2 for T states, up to 


~ 0. We have performed simulations of the dominant back¬ 
ground sources contributing to the invariant mass spec¬ 
trum. The uncorrelated background was obtained using Pythia 
generator and dimuons from correlated background sources: 
cc, bb and Drell-Yan, were simulated using both HELAC-Onia 
and Pythia generators. The estimated background level allows 
for y/i/i, ilf(2S), T(15), T(25) and T(35) measurements in the 
dimuon decay channel with good signal to background ratios. 

These simulations set the stage for further ones including, 
on the one hand, the detection of photon from P wave or tJc 
decay or from the production of a 7/i/i -H y pair, whose studies 
at low transverse momentum can provide important insight on 
the gluon transverse dynamics [51-54] and, on the other hand, 
the large combinatorial background typical of p+A and A+A 
collisions in which the study of excited quarkonium at AF¬ 
TER @LHC energies is of paramount importance [1, 11]. We 
note that the Delphes [55] framework seems particularly well 
suited to account for the photon detectability in such prospec¬ 
tive studies. 

Along our investigations, we have also noted that main 
source of dimuons around the T{nS ) masses is from the Drell- 
Yan process (see Fig. 6 (right)). This gives us great confidence 
that the corresponding cross section can easily be extracted in 
this mass region in p + p collisions, a fortiori with a vertex 
detector allowing for tagging the heavy-flavour muons. We 
therefore consider that the single-spin asymetries for Drell- 
Yan pair production can indeed be extracted using a light po¬ 
larised target. Motivations for such studies are discussed in 
[13, 14, 56]. Quarkonium polarisation measurement are of 
course also possible given the large statistical samples. 

As regards the case of p + A collisions, we have had a first 
look at the charged particle multiplicities as a function of the 
laboratory pseudo-rapidity. We have found out that, for all 
the possible hxed target modes, p+Ph, Pb-i-H, Pb-pPb, these 
are smaller than the ones reached in the collider modes where 
the LHCb was used (p-i-Pb and Pb-Hp at 5 TeV). We therefore 
believe that a detector with similar characteristics as compared 
to LHCb can very well be used in the fixed-target mode 

In view of the above, we have evaluated the impact -and its 
uncertainty- on the nuclear modification of the gluon densi¬ 
ties on prompt and non-prompt 7/i/i and T in form of Rp+pb- 
We have found that the measurements at backward rapidities 
allows one to search for the gluon antishadowing, the gluon 
EMC effect and even the Fermi motion effect on the gluons 
with unheard of statistical precisions. The statistics are large 
enough to perform such measurement with the i/'(25) and 
probably also with T(25) and T(35) allowing for thorough 
investigations of formation time effect of the meson propagat¬ 
ing in the nuclear matter. Overall, our results confirm the great 
potential of AFTER @LHC for heavy-quark and quarkonium 
physics. 


^ Our observation is obviously supported by the preliminary analysis of the 
LHCb-SMOG data taken during the pilot run of p* beam (Pb beam) on a 
Neon gas target from 2012 (2013) at a c.m.s energy of = 87 GeV 

(54 GeV) [57]. 
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